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We study properties of space-like monopole trajectories in the Maximal Abelian 
gauge of quenched SU(2) QCD at the finite temperature. We concentrate on in¬ 
frared monopole clusters which are responsible for the confinement properties of the 
theory. We determine numerically the effective action of the monopoles projected 
onto the three-dimensional time-slice. Then we derive the length distributions of 
the monopole loops and fix their entropy. 
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I. INTRODUCTION 


The interest in the Abelian monopoles in the non-Abelian gange theories is motivated by 
a central role of these objects in the dnal superconductor mechanism [3| of color confinement. 
The Abelian monopoles can be considered as particular conhgurations of gluon helds with 
magnetic quantum numbers. In pure non-Abelian gauge theories the Abelian monopoles 
do not exist at the classical level. However, these topological defects can successfully be 
identihed given a dynamical conhguration of the gluon helds in a particular Abelian gauge . 
There are maw Abelian gauges among which the most popular one is the Maximal Abelian 
(MA) gauge [3|. In this gauge the off-diagonal gluon helds are suppressed and short-ranged 
contrary to the diagonal (Abelian) gluon helds |^. There are many numerical experiments 
conhrming that Abelian degrees of freedom are responsible for the conhnement of color (for 
a review, see Ref. P). In particular, it was observed in Refs, 0 that the tension of 
the chromoelectric string is dominated by the Abelian monopole contributions. Moreover, 
the monopole condensate - which guarantees the formation of the chromoelectric string 
between the quarks - exists in the conhnement phase and disappears in the deconhnement 
phase [ 1 , 1 ^ . 

The trajectories of the Abelian monopoles form two diherent types of clusters. A typ ical 
conhguration contains a lot of hnite-sized clusters and one large percolating cluster |l0|, [ll| . 
The percolating cluster (or, the infrared (IR) cluster) occupies the whole lattice while the 
sizes of the other clusters have an ultraviolet nature. The monopole condensate corresponds 
to the so-called percolating (infrared) cluster of the monopole trajectories. The tension 
of the conhning string gets a dominant contribution from the IR monopole cluster 
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while the hnite-sized ultraviolet (UV) clusters do not play any role in the conhnement. 
Various properties of the UV and IR monopole clusters were investigated previously in 
Refs. 

At high temperatures the Abelian monopoles become static. In the high temperature 


phase the IR monopole cluster disappears ua lll| | and, consequently, the conhnement of 
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the static quarks is lost. Since the static currents do not play any role in conhnement we 
concentrate below on the spatial components of the IR monopole cluster. We investigate 
the action, the length distribution and the entropy of spatial components of the infrared 
monopole clusters. We follow Ref. where energy and entropy of the monopole currents 
were studied at zero temperature. Our preliminary results were reported in Ref. H- 

The plan of the paper is as follows. In Section E] we describe the model and provide 
the description of the monopole currents. The details of numerical simulations are also 
given in this Section. Section IIIII is devoted to the investigation of the Abelian monopole 
action obtained by the inverse Monte-Carlo method for the clusters of the spatially projected 
Abelian monopoles. In Section HVl we study the length distributions of the infrared clusters 
of the spatially projected monopole clusters. The knowledge of the monopole action and 
cluster distribution allows us to calculate the entropy of the spatial monopole currents which 
is discussed in Section El Our conclusions are presented in the last Section. 


II. MODEL 


We study pure SU(2) QCD with the standard Wilson lattice action for gluon helds, 

S((7) = -|H'^(7P, (1) 

P 

where {3 is the coupling constant, the sum goes over all plaquettes of the lattice, and Up = 
u,.,, = is the SU(2) plaquette constructed from link helds, We 

work in the MA gauge [3| dehned by the maximization of the lattice functional 

R = ^TT(a3U{s,^)a3U\s,^)^ , ( 2 ) 

S,/i 


with respect to the gauge transformations U{s,^) —> = Q{s)U{s, {s + fi). In 

the continuum limit the local condition of maximization (j2I) can be written in terms of the 
differential equation, — zA^) = 0. Both this condition and the functional (j21) 

are invariant under residual U(l) gauge transformations, = diag(e*^^®\ 

After the Abelian gauge is hxed we perform the projection of the non-Abelian gauge 
helds, onto the Abelian ones, Us,^'- 


U{s,g) 


(1 - |c(s,/i)n^/2 -c*(s,/i) \(u{s,g) 0 \ 

c(s,p) (1 - |c(s,/i)nR2 J 0 M*(s,/i)y’ 


( 3 ) 


where c(s, /i) corresponds to the charged (oh-diagonal) matter helds. 

As we have discussed above, the dominant information about the conhnement properties 
of the theory is located in the monopole conhgurations which are identihed with the help of 
the Abelian phases of the diagonal helds, 6s^^- The Abelian held strength 6^y{s) G (—dvr, dzr) 
is dehned on the lattice plaquettes by a link angle 0{s,fi) G [—vr, vr) as = 9{s,g.) -|- 

9{s + — 9{s -t- z>, p) — 9{s, v). The held strength 9^^[s) can be decomposed into two 

parts. 


9^,v{s) = 9^^{s) 27rm^^(s), 


( 4 ) 
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where 9^y{s 




G [—vr, vr) is interpreted as the electromagnetic flux through the plaquette and 
can be regarded as a number of the Dirac strings piercing the plaquette. 

The elementary monopole current can conventionally be constructed using the DeGrand- 
Toussaint 3 definition: 


2^^^upcrdu'^pcr(^^ T h)) (5) 

where d is the forward lattice derivative. The monopole current is defined on a link of the 
dual lattice and takes the values 0, ±1, ±2. Moreover the monopole current satisfies the 
conservation law automatically, 


= 0 , ( 6 ) 

where d' is the backward derivative on the dual lattice. 

The monopole current corresponds to the monopole charge defined on the scale of the 
elementary lattice spacing, a. Obviously, the scale a becomes smaller as we approach the 
continuum limit. In order to study the properties of the monopoles at fixed physical scales 
we use the so-called extended monopoles [l0|. The extended monopole is defined on a 
coarse sublattice with the lattice spacing b = na. Thus the construction of the extended 
monopoles corresponds to a block-spin transformation of the monopole currents with the 
scale factor n, 


n—1 

kfj,{ns + (n - l)fi +iu + jp + la). (7) 

i,j,l=0 

Since the time-like monopole currents are not essential for the confinement properties 
we concentrate on the spatial components of the currents. Namely, we investigate spatially 
projected currents. 


S4=0 


1,2,3, 


( 8 ) 


which are integer-valued and closed. 

Technically, we generate 2000-10000 configurations of the SU(2) gauge held, U, for (3 = 
2.3 ~ 2.6 on the lattices x Lt, with Lg = 24, 32,48, 72 and = 4, 6, 8,12,16. The number 
of generated configuration depends on the value of f3 and lattice volume. We fix the gauge 
with the help of the usual iterative algorithm. In this paper we used the same methods 
as in the zero-temperature case studied in Ref. Q. Thus we refer an interested reader to 
Ref. 0 for a more detailed description of the numerical procedures. Below we concentrate 
on the description of the numerical results. 


III. MONOPOLE ACTION 

In what follows we discuss an effective model of the monopole currents corresponding to 
pure SU(2) QCD. Formally, we get this effective model through the gauge fixing procedure 
applied to the original model. Then we integrate out all degrees of freedom but the monopole 
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ones. An effective monopole action is related to the original non-Abelian action S'[t/] as 
follows: 

/ CXD 

OTi(A')Arp(C/)e-“I^I= (n E ) ■ (9) 

kfj,{s) = — 00 s 

We omit irrelevant constant terms in front of the partition function. The term 6{X) repre¬ 
sents the gauge-hxing condition and Afp{U) is the corresponding Faddeev-Popov determi¬ 
nant. As we have discussed above, the MA gauge hxing condition is given by a maximization 
of the functional (j2|) and therefore the local condition X = 0, implied in Eq. (jH)), is used 
here as a formal simplihed notation. 

Numerically, the monopole action of the 3D projected IR monopole clusters can be dehned 
using the inverse Monte-Carlo method . The action is represented in a truncated form |^, 
E3 as a sum of the m-point (m > 2) operators Sp 

Smon[K] = y^ f^Si[K], ( 10 ) 

where fi are the coupling constants. Following Ref. |l^ we adopt only the two-point inter¬ 
actions in the monopole action. Si Kds)K,(s'). 

Similarly to the AD case we hnd that the monopole action of the spatially projected 
currents is proportional with a good accuracy to the length L[K] of the monopole loop K, 

5'mon[A'] ^ foL[K] + const. (11) 

The important property of the monopole action is that the couplings fi are the functions of 
the scale b = na, Eq. o, at which the monopole charge is dehned. To illustrate this fact 
we show the dependence of the coupling constant /q on 6 = na{j3) in Figure [TJ 
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FIG. 1: The coefficient /o of the monopole action m vs. the scale parameter b for the lattice 
sizes Lg X 6, Ls = 48, 72 and blocking factors, re = 1... 9, at temperatures (a) T = 0.8 Tc and (b) 
T = 0.96 Tc. 


From Figure Q one observes the almost perfect scaling: the parameter b does not de¬ 
pend on the parameters n and a separately. The action is near to the renormalized tra¬ 
jectory which corresponds to the continuum effective action. Moreover, the result does 
not depend on the spatial extension of the lattice, Lg. Thus the action of the spatially 
projected monopole current shows the scaling similarly to the action of the unprojected 
monopoles [1, |lg . 
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IV. LENGTH DISTRIBUTION 

The length distribution of the spatially-projected monopole clusters is shown in Figure|21 
In the confinement phase, T < Tc, only the infrared part of the distribution is shown. One 




FIG. 2: The distributions of the spatial monopole currents at various temperatures. For the low 
temperatures, T < Tc, the UV-part of the distributions is not shown. 

can see that the length (in physical units) of the monopole trajectory belonging to the 
percolating cluster becomes shorter as the temperature increases. This fact is expected 
because the monopole condensate is ’’evaporating” as temperature increases towards the 
transition point, and, therefore, the infrared part of the monopole currents should be more 
and more diluted. 

At T > Tc the percolating cluster of the spatially projected currents disappears, and, 
consequently, the confinement of quarks is lost. The behavior of the elementary and blocked 
currents is qualitatively the same. According to Ref. Q the length of the AD IR monopole 
currents in the hnite volume V is distributed with the Gaussian law, which is in the finite- 
temperature case can be formulated as follows: 

D^^{L) oc exp{-a{b, V)L‘^ + 7 ( 6 , T)L} . (12) 

The length distribution function, D{L), is proportional to the weight with which the partic¬ 
ular trajectory of the length L contributes to the partition function. In Eq. m we neglect 
a power-law prefactor, l/T’’" with r ~ 3, which is essential for the distribution of the infrared 
clusters. 

The Gaussian form of the distribution (HI means that the clusters have the typical length 

Lmax ='l{h,T)/2a{h,V), (13) 

where V is the three-dimensional volume. The coefficient a plays a role of the infrared 
cut-off which emerges due to the hnite volume. In other words, the length of the monopole 
trajectory in an infrared cluster is restricted by the lattice boundary. However, since the 
cluster is infrared the length of the monopole trajectory in this cluster must be proportional 
to the total volume, L^ax oc V. The linear part of the distribution m gets contribution 
from the monopole action and the monopole entropy (we discuss this issue below). Therefore 
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the coefficient 7 should not depend on the volume in the thermodynamic limit. Thus, we 
expect 


a{h,V) = A{h)lV, 


(14) 


where A{h) is a certain function of the scale parameter b. One may suggest that the pa¬ 
rameter A should not signihcantly depend on the temperature T since the factor is more 
kinematical than dynamical. The temperature influences the dynamical characteristics of 
the monopoles such as the effective three-dimensional action. The effective monopole ac¬ 
tion contributes to the coefficient 7 and, as a consequence, the temperature influences the 
projected monopole density via the 7 -coefficient. 

Using Eqs. (USUI one can obtain that the monopole density in the infrared cluster is 
hnite in the thermodynamic limit and is given by the formula 

Ljnax _ 'yiPi 

PIR = y = 2 71(5) • 



We £t the numerically obtained distributions of the 3D projected currents by the func¬ 
tion m and then use a bootstrap method^ to estimate the statistical errors of the htting 
parameters. In Figure El we show the coupling constant 7 ( 6 , T) as a function of the scale 
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FIG. 3: The same as in Figure Q but for the coefficient 7 of the distribution m- 


parameter, b, at temperatures T = 0.8 Tc and T = 0.96 Tc. Again, as in the case of the pa¬ 
rameter /o. Figured we observe the volume independence and the fe-scaling of the results. 
In a small 6 -region we find that 7 oc 6 ^ with 7 ~ 3 for low temperatures, T ~ O.STc, whereas 
7 ~ 2 for T —>■ Tf.. The data show a good 6 -scahng and also is independent of the volume 
similarly to the monopole action. 

The numerical values of the parameter A are shown in Figure d The parameter A is 
independent of the lattice volume, indicating that in the thermodynamic limit the coefficient 
a in the Gaussian distribution m vanishes. 


^ A detailed description of the corresponding bootstrap method is given in Ref. Ih 
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FIG. 4: The same as in Figure ^ but for the ratio A(b), Eq. (1141) . 
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FIG. 5; The same as in Figure ^ but for the monopole density p corresponding to the infrared 
cluster, Eq. m- 


Using Eq. m we calculate the monopole density corresponding to the infrared cluster. 
The density is shown in Figure |S1 The density diminishes as the scale factor b increases, 
while at small b the density shows a plateau. As the temperature increases the density (at 
a hxed value of b) becomes smaller. 

Note that the conhning non-Abelian objects have a finite size (in physical units). These 
objects are identihed as the Abelian monopoles in the Abelian gauge. The monopoles are 
detected using the Gauss theorem applied to the magnetic held coming outside the cube of 
the size b^. The conhning non-Abelian objects have a typical size which is associated with 
the size of the monopole core Q , rmon ~ 0.05 fm. If 6 < rmon then the monopole cube is too 
small to detect the charge of much larger monopole and the monopole density - measured 
using the Gauss law - is vanishingly small. Indeed, one can see that the monopole density 
in Figure El has a tendency to diminish at smaller by/a <0.1. 

At the critical temperature, T = Tc, the AD IR monopole cluster disappears and we 
expect a similar behavior for the 3D projected IR cluster. This implies, that the parameter 
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FIG. 6 : The coefficient 7 vs. temperature T for various temporal extensions of the lattice and for 
various blacking factors n. 
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FIG. 7: Examples of the fits of the 7 parameter as the function of the temperature. 


7 ( 6 , T) must become a (non-local) order parameter, it must vanish at the critical point. 
One can reach this conclusion by noticing that the parameter 7 ( 6 , T) is proportional to the 
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n 

48^ X 6 

48^ X 8 

72^ X 8 

1 

0.64(15) 

0.76(2) 

- 

2 

0.62(8) 

0.70(16) 

0.48(3) 

3 

0.34(6) 

0.55(7) 

0.30(2) 

4 

0 . 22 ( 2 ) 

0.36(6) 

0.18(3) 

6 

0 . 11 ( 2 ) 

0 . 20 ( 2 ) 

- 


TABLE I: The ’’critical exponents” 5 - obtained with the help of the fit m - for various lattices 

X Lt and extensions n. 

monopole density ()15jl (which vanishes at T = Tc), and that the factor A is unlikely to be 
divergent at the critical temperature (what can also be deduced from Figures |3)). 

We show that the quantity 7 is indeed an order parameter in Figure IHKa) which depicts 
7 for elementary (n = 1 ) monopoles as a function of temperature for various temporal 
extensions of the lattice. The behavior of the 7 -parameter in the vicinity of the phase 
transition point depends on the value of the temporal extension Lt. However, the value of 
the 7 -parameter at the critical temperature, 7 (T —> T^) —>■ 0, is universal with respect to 
Lt- Moreover, one can observe in Figures ini(b,c,d) - which correspond to the extensions 
n = 2,3, 4, respectively - that the 7 -coefhcient for the extended monopoles is also vanishing 
at T = Tc- 

To characterize the critical behavior of the parameter 7 in the vicinity of the phase 
transition we have htted this parameter by the function 

7“(6,r) = C,. (l-L)'. T<n, (16) 

where 6 and C.y are the htting parameters. We performed the hts for various lattices x Lt 
and extensions n. The results for the ’’critical exponent” 6 are shown in Table HI From this 
table one notices that the quantity 6 is not universal: it depends not only on the extension of 
the monopole blocking but also it also depends on the lattice volume. Moreover, the larger 
extension n the steeper behavior of 7 in the vicinity of the phase transition is. 

One should add a word of caution here. The £t results shown in Table H] are crucially 
dependent of the T/Tc = 0.98,0.99 points (as one can see from Figures [7Ka),(b)), which 
are very close to the phase transition. Since the transition is of the second order then the 
hnite-volume effects must be strong and the results of hts may quantitatively be incorrect 
(although the results presented in Figures [7Ka),(b) must qualitatively be correct). 


V. MONOPOLE ENTROPY 

Apart from the hnite-volume ehect, the distribution m has contributions from the 
energy and the entropy. As seen above, the action contribution is proportional to 
The entropy contribution is proportional to (with /i > 0) for sufficiently large monopole 
lengths, L. Thus, the entropy factor, /r, is 


A = exp{/i, + 7 } . 


(17) 
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We determine the entropy using Eq. (ED- The numerical results for the entropy factor 
/i(6, T) are shown in Figures|Hlfor various temperatures, lattice volumes and blocking factors. 
One can see that the entropy factor /x scales as the function of b, as expected. 

In order to understand the meaning of the data shown in Figure |H] we note that if the 
monopoles are randomly walking on a 3D hypercubic lattice then we should get a dehnite 
value for the entropy factor, /i = 5. This is because at each site there exist hve choices 
for the monopole current to go further. One can see that far from the phase transition, 
T < 0.96Tc the entropy factor /x indeed tends to the /x = 5 plateau at moderately small 
values of 6 ~ 0.4... 1. At yet smaller b the entropy gets bigger than random walk value 
/i > 5 because in this region the inverse Monte-Carlo method with the truncated quadratic 
monopole action does not work well jl^. Thus, the value of the constant /o ^ dehned in 
Eq. (ITT|l - can not be obtained correctly. 

At large b the entropy factor drops down with the increase of the factor b. This feature 
is independent of the temperature. In the zero temperature case jl^ the entropy factor /x 
approaches unity in the 6 —> cx) limit. This feature is difficult to observe from our data since 
the information about the entropy at large values of the blocking size b is not available. 


VI. CONCLUSION 

The distributions of the spatially-projected infrared monopole currents of various block¬ 
ing sizes, n were studied on the lattices with different spacings, a, and volumes, x Lf. 
We find that the distributions can be described by a gaussian anzatz with a good accuracy. 
The anzatz contains two important terms: (i) the linear term, which contains information 
about the energy and entropy of the monopole currents; and (ii) the quadratic term, which 
appears due to finite-volume and which suppresses large infrared clusters. The linear term 
is independent of the lattice volume while the quadratic term is inversely proportional to the 
volume. Moreover, the linear term is a (non-local) order parameter for the deconfinement 
phase transition. 

To get the entropy of the spatially-projected currents we studied the action of the 
monopoles belonging to the infrared monopole clusters of the spatially-projected currents 
using an inverse Monte-Carlo method. We show that the entropy factor has a plateau at 
sufficiently small values of b and at T < 0.96Tc. A reason for the temperature restriction of 
our result is that our analysis may not be valid close to the second order phase transition 
point because of the increase of correlation lengthes at (and, consequently, because of strong 
finite-volume effects) T ^ T^. At 6 > 1 the entropy is a descending function oi b = na, 
indicating that the effective degrees of freedom of the projected and blocked monopoles 
are getting smaller as the blocking scale b increases. This effect is very similar to the zero 
temperature case, in which the monopole motion corresponds to the classical picture: the 
monopole with the large blocking size b becomes a macroscopic object and the motion of 
such a monopole gets close to a straight line. 


11 



7 

5 

4 

3 

2 

1 

0 



• 



Random walk in D=4 : 


T=0.8T 


• 


C 



-A - 



A 

□ 

Random walk in D=3 : 


♦ 

▼ fi 




♦ 





s 

7^x6 : 

□ 

n=2,3,4,6,9 

^ □ : 

: 4^x6 : 

• 

n=2 A n=3 



▼ 

n=4 ♦ n=6 

^ n=8 


0 




(b) 


7 

/^6 

5 

4 

3 

2 

1 

0 


: i .: 


Random walk in D=4 : 


T=0.9T 

C 

1 Random walk in D=3 : 

: 

i 

i 


- 72^x6: □ 

□ 

n=2,3,4,6^ - 

: 4^x6: • 

n=2 A n=3 

▼ 

n=4 ♦ n=6 ^ n=8 


0 


b<5 


(c) 


1/2 3 


7 

/^6 

5 

4 

3 

2 

1 

0 


Random walk in D=4 


T=0.96T 


Random walk in D=3 




1 □ 

7^x6\ □ n=2,3,4,S,9 
4^x6: • n=2 A n=3 

T n=4 ♦ n=6 * n=8 


1 




1/2 


(d) 


7 

5 

4 

3 

2 

1 

0 , 


Random walk in D=4 


T=0.98T 


Random walk in D=3 


72^x6: □ n=2,3,4,6,9 
4^x6: • n=2 A n=3 

T n=4 ♦ n=6 ^ n=8 


1 




1/2 3 


(e) 


7 

/^6 

5 

4 

3 

2 

1 

0 


Random walk in D=4 

T=0.99T 


□ 

"'ii* I 

72^x6: □ n=2,3,4,S,9 
4^x6: • n=2 A n=3 

T n=4 ♦ n=6 


Random walk in D=3 

1 □ 


* n=8 


1 




1/2 


(f) 


FIG. 8; Entropy factor of the spatially projected monopole currents as the function of the scale b 
at various temperatures. 
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